Does the mass of a black hole decrease due to the accretion of phantom energy? 
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According to Babichev et al, the accretion of a phantom test fluid onto a Schwarzschild black hole 
will induce the mass of the black hole to decrease, however the backreaction was ignored in their 
calculation. Using new exact solutions describing black holes in a background Friedmann- Robertson- 
Walker universe, we find that the physical black hole mass may instead increase due to the accretion 
of phantom energy. If this is the case, and the future universe is dominated by phantom dark energy, 
the black hole apparent horizon and the cosmic apparent horizon will eventually coincide and, after 
that, the black hole singularity will become naked in finite comoving time before the Big Rip occurs, 
violating the Cosmic Censorship Conjecture. 



I. INTRODUCTION 



The first year WMAP data combined with the 2dF 
galaxy survey and the supernova la data favor the phan- 
tom energy equation of state of the cosmic fluid w = 
P/p < — 1 (where p and P are the cosmic energy den- 
sity and pressure, respectively) over the Einstein cos- 
mological constant [w — —1) and the quintessence field 
(w > — 1) 1, 2]. The interpretation of this result brings 
about the possibility that dark energy could be in the 
form of a phantom field [!, 0, [f| . Compared to the La- 
grangian density of a quintessence field ^quintessence = 
Tjduifid^cf) + V (<j)), the Lagrangian density of the phan- 
tom field has a kinetic energy term with the "wrong" 
sign, i.e., ^phantom = -|9 M </>9 M + V (</>). This"wrong" 
sign endows the phantom model of dark energy with two 
surprising properties. First, the density of phantom mat- 
ter increases with the expansion of the universe in such a 
way that both the scale factor and the phantom energy 
density can become infinite in a finite comoving time. 
If this happens, every object in the universe, including 
galaxies, stars, the earth and nuclei will be torn apart in 
a finite future in a so-called Big Rip singularity [1, B 13 • 
A second surprising property of the phantom field was 
found by Babichev, Dokuchaev and Eroshenko Q, who 
showed that accretion of phantom energy will induce the 
mass of black holes to decrease. They concluded that the 
masses of all black holes will vanish before the Big Rip 
is reached. However, in their calculations, they ignored 
the backreaction of phantom matter on the black hole 
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metric. In a low matter density background, this effect 
can be safely ignored, but if the matter density of the 
background is large, for example if it becomes compara- 
ble to the black hole density ~ 1/Mg (Mq is the black 
hole mass), the metric describing this black hole will be 
modified significantly; hence, the backreaction must be 
accounted for in this case. 

The purpose of the present paper is to point out that, 
taking into account the backreaction, the physical black 
hole mass may not decrease but will rather increase due 
to the accretion of phantom matter. The size of the black 
hole apparent horizon (AH) increases, while that of the 
cosmic AH decreases, during the evolution of the universe 
to a Big Rip. There exists a moment t* at which the 
black hole AH coincides with the cosmological AH, and 
after which the black hole singularity becomes naked. A 
naked singularity is forbidden by the Cosmic Censorship 
Conjecture [9(; in a broad sense, if the latter is correct, 
it seems that the Cosmic Censor disallows the existence 
in nature of not only naked singularities, but also of the 
cryptic phantom field. 

In order to study this issue we use new exact solu- 
tions of the Einstein field equations describing dynam- 
ical black holes embedded in a Friedmann-Robertson- 
Walker (FRW) universe driven by phantom energy and 
accreting this cosmic fluid. Independent motivation for 
the development of dynamical black hole solutions comes 
from various areas of theoretical physics including black 
hole thermodynamics, quantum gravity, and mathemat- 
ical physics [IE m E H Q El ■ Similar issues about 
the accretion of phantom energy onto wormholcs were 
discussed in [lj|. 

We use units in which c = G = h = 1 and this paper 
is organized as follows: in Sec. II we discuss how the 
McVittie solution can not describe a black hole embedded 
in a FRW background universe. Then, we turn to new 
solutions and show that a new metric found (but not 
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studied in detail) in [l7| can describe a dynamical black 
hole in a universe driven by dark matter, radiation, or 
dark energy. In Sec. Ill, using the solution of [13], we 
derive the result that the black hole mass will increase 
due to the accretion of phantom matter. Sees. IV and VI 
find the same result for qualitatively different classes of 
exact solutions, while Sec. VI contains a discussion and 
the conclusions. 



II. NEW COSMOLOGICAL BLACK HOLE 
SOLUTIONS 

To investigate the relation between cosmic expansion 
and local physics, McVittie introduced in 1933 his well- 
known solution \li 
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where dfl 2 is the line element on the unit 2-sphere and 
the constant Mq denotes the physical mass of the black 
hole when a(t) =const. (this can be seen by rescaling the 
radial coordinate as f = ar). The scale factor a(t) is an 
arbitrary (positive) function of the comoving time t. This 
metric reduces to the familiar spatially flat FRW metric 
as M — > 0. So at a first glance, the metric seems to repre- 
sent a Schwarzschild black hole embedded in a spatially 
flat FRW universe [lj|. However it is argued that, in gen- 
eral, the McVittie solution can not describe a black hole 
because it is singular on the 2-sphere r = Mq/2 (which re- 
duces to the Schwarzschild horizon if a =const.) [2fJ and 
this singularity is spacelike [2lj |. It is also claimed that 
the McVittie metric may describe a point mass located 
at r = and embedded in a FRW universe. However, 
this point mass is, in general, surrounded by the singu- 
larity at r — Mo/2, which is difficult to interpret. This 
singularity was studied in Refs. (20|. Nolan [22| showed 
that this is a weak singularity in the sense that the vol- 
ume of an object approaching the r = Mo/2 surface is 
not shrunk to zero, and therefore the energy density of 
the cosmic fluid remains finite. On the other hand, the 
pressure diverges at r = Mo/ 2 together with the Ricci 
scalar R [20l . [23j]. So, while the McVittie metric rep- 
resents some kind of strongly gravitating central object, 
its physical interpretation is not totally clear and remains 
the subject of debate [2(| EH- The only exception is the 
Schwarzschild- de Sitter (or Kottler) metric, which cor- 
responds to a special case of the McVittie metric and 
represents a genuine black hole of constant mass embed- 
ded in a de Sitter background [HEilEl- 

Because the McVittie solution can not describe a black 
hole in a FRW universe (with the exception of the 
Schwarzschild-de Sitter one), we turn our attention to 



the new solutions of Ref. [l7| , described by the general- 
ization of the McVittie metric 1241 
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in the background of an imperfect fluid with a radial heat 
flux and, possibly, a radial mass flow simulating accretion 
onto a black hole embedded in a generic FRW universe. 

In the McVittie metric, the mass coefficient M is con- 
stant, which expresses the assumption (hereafter referred 
to as the "McVittie condition") that G\ — 0, where G^ v 
is the Einstein tensor, and therefore to vanishing radial 
mass flow onto the central object, Tq 1 = 0, and to no ac- 
cretion. M(i) corresponds [l?]] to the Hawking-Hayward 
quasi-local mass of the metric (Eq. [2]) and, therefore, the 
McVittie no-accretion condition corresponds to constant 
physical mass of the central object, which is static even 
when the universe expands or contracts. Here instead, 
following [13] , the black hole is allowed to accrete cosmic 
fluid and M(t) changes with time. 

Two distinct processes concur in the time variation of 
the black hole mass: the first is the radial flow of energy 
onto the black hole due to the radial motion of the cos- 
mic fluid toward the horizon (this is a simplified version 
of the accretion process considered by Babichev et al.), 
possibly supplemented by a radial heat flux, which neces- 
sarily requires an imperfect fluid description. The second 
process is the variation of the mass enclosed by the black 
hole apparent horizon due to the motion of the horizon 
itself; this process is absent in the McVittie solution (in 
which the singular surface r = M /2 does not expand) 
and in the static black hole considered by Babichev et 
al, but is present for any dynamical horizon even in the 
absence of the first accretion process. However, unless 
the total mass of the central object is zero (which could 
happen, e.g., for a spherical shell of exotic matter), a ra- 
dial energy flow is to be expected due to its gravitational 
pull. The solutions of [13] exhibit accretion onto the cen- 
tral black hole due to radial mass flow and/or to a radial 
heat flux. 

The components of the Einstein tensor for the met- 
ric © are [l3] 
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where A = 1 + B = 1 - and m{t) = M(t)/a(t). 
Assuming that there is no mass flow but that there is a 
radial heat flux in the stress-energy tensor 

T^v = {P + p) u^u v + Pg^ v + q^u v + q u u^ , (6) 

where = (|-goo|~ 1/2 . °> °> °) and ^ = (0,9,0,0), the 
Einstein equations are 
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For the constant equation of state P = w p (w =const.), 
combined with the (0,0) and (1,1) Einstein equations to 
obtain G\ + w Gq — 0, one derives 
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is non-negative. 

Thanks to the particular form of the metric coefficient 
M (t) , this solution is conformal to a Schwarzschild black 
hole with conformal factor equal to the scale factor of 
the background universe a(t) [l a . \1 7\, a feature shared 
with the Sultana-Dyer solution 25] (c/. also the solu- 
tions of Ref. [23j]). The mass coefficient M(t) coincides 
with the Hawking-Hayward quasi-local mass of the space- 
time (EqJ2), which should be identified with the physical 
black hole mass 15]. The fact that the latter increases in 
an expanding universe describes the constant accretion 
due to the radial energy flux from the surrounding fluid. 
The McVittie condition G\ = equivalent to M = that 
forbids accretion onto the central object in the McVittie 
metric is clearly not imposed in the solutions of [13] that 
we use here, allowing accretion to occur. 

These solutions describe a cosmological black hole em- 
bedded in a FRW background and present advantages 
over the previous Sultana-Dyer [25| and McClure-Dyer 
[23[ solutions, in the sense that both the energy density 
and the pressure are finite near the black hole horizon, 
and the energy density is positive-definite [l7l ]. 

More general solutions with a different form of the 
mass coefficient and no conformal Killing horizon are con- 
sidered later in Sec. IV. 
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which expresses the fact that an ingoing flow (q < 0) 
corresponds to an increase of the quasi-local mass M > 0, 
and vice-versa. 

As r — > +00, the scale factor of the universe is given 
by the familiar expressions 



a(t) = a (t np -t) 3 ^+ 1 i if w < -1 , (12) 
a(t) = a e Hot if w = -1 , (13) 



a(t) = a i 3 <™ +1) if w > -1 



(14) 



where t r i v denotes the Big Rip. It is easy to see that 
the previous form of a(t) with the choice M(t) = Mo a(t) 
(with M a positive constant) and the corresponding heat 
flux density 
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constitutes a solution of the full Einstein equations. The 
energy density and pressure depend on both the radius 
and the comoving time, p = p (t,r), P = P (t,r), and 
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III. VIOLATION OF COSMIC CENSORSHIP 

The exact solution describing a black hole embedded in 
a FRW universe that we adopt in this section corresponds 
to the choice M (t) = M a(t) and a(t) given by Eq. (Tl2|) . 
and it can be written as 
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in isotropic coordinates. If a = 1, then r = Mo/2 is the 
event horizon of the Schwarzschild black hole; however, 
if a is not constant, r = Mo/2 is no longer the event 
horizon but merely a 2-spherc enclosed within the AH 
and a conformal Killing horizon. 

In order to investigate the evolution of the black hole 
mass and horizon, it is convenient to rewrite Eq. (|17p 

using the areal radius r = r(l + ^) . Eq. (fT7|) then 
becomes 



l-^W + a'Tl-^l dr 2 



ds' = - 

+a 2 f 2 dn 2 . (18) 
Further setting R = af, Eq. (jT5J) is rewritten in the 
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Painleve-Gullstrand form as 
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The coordinate system (i, R, 6, <fi) is not orthogonal; 
the cross-term dtdR can be eliminated by introducing 
the new time coordinate t defined by 



dt 



1 



dt 



HR 



-dR 



(20) 



where F (t(t, R), R) is an integrating factor which always 
exists and solves the partial differential equation 
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in order to make dt an exact differential. The metric {19 
is cast in the Nolan gauge form 
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where H (t) and F(t, R) are now functions of i and R. 
The Schwarzschild solution in Schwarzschild coordinates 
is recovered by setting a = 1 (and, consequently, H = 
and R — f), while the Schwarzschild-de Sitter solution 
is recovered by using Eq. (1) and Eq. (13); further, the 
latter reduces to the de Sitter solution (in static coordi- 
nates) if one also sets M = 0. 

In the Schwarzschild-de Sitter spacetime, there exists 
a black hole apparent horizon (AH) and a cosmic AH. 
The metric (|22p describes a spherically symmetric and 
dynamical black hole in a FRW background more gen- 
eral than de Sitter. In this case, the event horizon may 
not be well-defined, but the AH always exists. The AH 
is a marginally trapped surface with vanishing expan- 
sion and has been argued to be a causal horizon for a 
dynamical spacetime. The AH is associated with the 
Hawking temperature, gravitational entropy and other 
thermodynamical aspects [l^, HH, EH EH • The first law 
of thermodynamics for the AH has been derived not only 
in general relativity but also in several other theories of 
gravity including the Lovelock, nonlinear, scalar-tensor, 
and braneworld theories [H, [23, [H, [H, S [H H| . In 



view of this point, in order to investigate the evolution 
of the black hole mass, we calculate the energy contained 
inside the AH. To this end, let us first calculate the radius 
of the black hole AH: Eq. $2$ yields 
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as the equation of the AH (this expression appears also in 
Ref. [l4| , and it can also be derived from Eq. (fTT|) by re- 
placing r with R) . We discard the solution corresponding 
to the lower sign and to a negative radius, then Eq. 
has the two physical roots 
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where an overdot denotes differentiation with respect to 
the comoving time t and the explicit time-dependence of 
H and d contains the time evolution of the AH in which 
we are interested. R c and Rb identify the cosmic AH 
and the black hole AH, respectively. When the radius 
R crosses the black hole AH R},, the goo and gn terms 
in the line element (f2"2"|) change sign, revealing that the 
AH is a one-way membrane for particles falling into the 
black hole. A calculation of the null surface of equation 

9^ 6^ dl» = shows that the AH is a nul1 S-sphere- The 
coordinate transformations yield the correspondence 

r = M /2 f = 2M R = 2M a = 2m H {t)(26) 

where ttih — M is the Hawking-Hayward quasi-local 
mass (see the end of this section for a discussion). So, 
the conformal Killing horizon r = M /2 is no longer the 
black hole event horizon but is instead a 2-sphere con- 
tained inside the AH. In fact, it follows from Eq. (j2"3")l 
(with the upper sign) that the roots of this equation sat- 
isfy R Cy b > 2M a — 2itih- Moreover, the energy density 
p and the pressure P are both regular on the AH and p 
is positive-definite. 

We are interested in the accretion of phantom energy 
onto the black hole; the black hole AH and the cosmic 
AH computed from Eq. (fT2")) are plotted in Fig. ^ for the 
parameter choice eto = 1, Mq = 1, t r i p — 1, and w = —2. 
As is clear from this figure, in our solution the black 
hole AH increases while the cosmic AH decreases as the 
universe evolves. If the future universe is dominated by 
phantom energy, then there is a time t* at which these 
horizons coincide, which is given by d(i*) = -gjj^, or 
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At the comoving time , the critical AH radius is R cr = 
^jj. At times t > i* both horizons disappear and the 
black hole singularity becomes a naked singularity. We 
stress that this happens before the Big Rip is reached. 
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The appropriate notion of mass for a cosmological 
black hole is the Hawking-Hayward quasi-local mass 
which, for the metric ([2| reduces to TOjj(f) = M(t) and, 
for the cosmological black hole (fTT)) , to 



m H (t)=M a(t) 



(28) 



For our solution this is always increasing in an expand- 
ing universe, even when the black hole AH disappears 
(the computation of this mass does not depend on the 
presence of an AH), a fact that can be interpreted as the 
energy flow continuing to fall onto the naked singularity 
after the horizon has disappeared. 

The conclusion that the black hole AH and the Hawk- 
ing mass increase in a phantom-dominated universe run- 
ning toward the Big Rip is the opposite of what Babichev 
et al. 8] find by extrapolating rigorous results ob- 
tained for the accretion of a phantom test fluid onto 
a static Schwarzschild black hole. They find that the 
Schwarzschild mass decreases because the accretion rate 
to oc (P + p) < 0. However, this result can not be gener- 
alized directly to a cosmological, time-dependent, black 
hole. The reason for the discrepancy is that the Hawking- 
Hayward quasi-local mass to# = M (which implies the 
subtraction of a background energy and reduces to the 
Schwarzschild mass when a = 1) depends only on the 
energy density p but not on the pressure P [33|]. Since 
p > 0, it does not matter how negative P is, and tojj 
is always positive. The Schwarzschild mass notion is not 
appropriate, nor defined, for a black hole in an asymp- 
totically FRW spacetime. 

The quasi-local mass has been proved to be superior 
to other notions of mass, which have generated some 
confusion in the literature, in the study of the Hawk- 
ing temperature of cosmological black holes [l5| . Never- 
theless, it may be interesting to try out a different no- 
tion of mass. An obvious candidate is the (generalized) 
Misner- Sharp mass that has been used in black hole ther- 
modynamics in relation with energy flows through the 
AH [13, [Hi- While we do not want to fuel a debate 
on which is the most convenient notion of mass, and for 
which purposes, we present the time evolution of this 
generalized Misner-Sharp mass contained inside the AH. 
For a spherically symmetric space-time with line element 
ds 2 = h^dx^dx" + R 2 dil 2 , the generalized Misner-Sharp 
mass is M M s =R(1- h» u R^Ru) /20|. At the AH, it 
is h^ v R 4l R. u — and the generalized Misner-Sharp mass 
inside the AH (also computed in Ref. [13]) is simply 
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for all times for which the expression (|25|) of Rb is 
well-defined. Both the quasi-local and the general- 
ized Misner-Sharp mass of the black hole increase with 
time: this is at odds with the decreasing mass prop- 
erty found in Ref. [1] (the difference between Misner- 
Sharp and Hawking-Hayward masses can be seen also in 
the Schwarzschild-de Sitter metric). Eqs. (|2"5|) and (j2"9")l 
show that a discussion of the time evolution of black hole 
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FIG. 1: The black hole apparent horizon (solid line) increases 
while the cosmic apparent horizon (dotted line) decreases with 
the accretion of phantom energy onto the black hole. There is 
a time at which the two horizons coincide and after which both 
horizons disappear and the black hole singularity becomes a 
naked singularity. The plots correspond to the equation of 
state parameter w = —2. 



masses is meaningless until the correct notion of mass 
is identified, and the choice is non-trivial for a time- 
dependent cosmological black hole. 



IV. A MORE GENERAL CLASS OF 
SOLUTIONS 



One may wonder whether the conclusions reached in 
the previous section are due to the peculiar form of the 
solution employed, in particular the choice of Hawking- 
Hayward quasi-local mass scaling as the scale factor, 
M(t) = Moa(t). The answer is negative, as we shall see 
below. The exact solution studied in the previous section 
can be generalized to an arbitrary form of the Hawking 
mass inn = M(t) [l7| (in this case, the 2-sphere r = M/2 
is no longer a Killing horizon). Consider again Eq. @ 
with a general form of M(t); by defining the areal radius 

f = r ^1 + 2ra(t) ) as before and m(f) = M(t)/a(t), the 
metric @ is rewritten as 



ds 2 
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Introducing the comoving R = af, in terms of which 
dr = - Hfdt, Eq. (30]) is turned into the Painleve- 



6 



Gullstrand form 
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mass and that of the scale factor of the cosmic substra- 
tum. The vanishing of this factor can be regarded as 
the analog of the condition for stationary accretion in a 
time-dependent background (keeping with the spirit of 
the quasi-local mass definition in which a background 
quantity is subtracted) . Then, the solution considered in 
the previous section corresponds to stationary accretion 
relative to the FRW background. 

By discarding the lower sign as usual, Eq. becomes 



'- | dtdR . 
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Further setting A(t, R) = 1 - 2M/R and C(t, R) 
and defining the time coordinate t by 
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where F (t, R) is again an integrating factor that makes 
dt an exact differential, one has 



(37) 

Due to the general form of the function M(t) (not scaling 
as M(t) = M a(t)), now the coefficient (l + ^) appears, 
in which r = r(R). Due to this dependence, Eq. (|3T|) can 
not be regarded as a quadratic algebraic equation for the 
AHs radii as before. Nevertheless, it is instructive to 
treat it formally as such, obtaining the formal roots 
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The cross terms in dR dt cancel out and we are left with 
the line element in the Nolan gauge 
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Since r depends on R, this is only an implicit equation for 
the cosmic and black hole AH radii R c ,b, but it suffices for 
our purposes. When the argument of the square root is 
positive, there are a cosmic AH at R c and a black hole AH 
at Rb while, if this argument becomes negative because 
the time derivative d increases without limit as the Big 
Rip is approached, these AHs coincide at the common 

critical radius \J 2 jf, and then disappear leaving behind 

a naked singularity. This situation occurs, e.g., for the 
simple choice of the mass function 



+R 2 dn 2 
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3\w + l\a«- 



t„ p and 



t) 3 <™ + 1 ' 



Since the Hawking mass is M(t) = m(t)a(t), the expres- 
sion on the left hand side can be written as 



(39) 



(40) 



3M-1 



HR + M 



1 777 \ 




( 1+ 27) 









(36) 



in which the factor M (l + measures the deviation 
of the radius from 2M (with r > m/2 ^ R > 2M 

M (1 + ^f) > 2M), while the factor - H^j is the 

difference between the percent variation of the black hole 



and the heat flux (fTTj) is finite on the AH. If a < 
(which is compatible with both requirements a > 1 and 
w < —1), the horizon disappears but the central mass 
does not, leaving a naked singularity. 



V. ANOTHER DYNAMICAL PHANTOM 
BLACK HOLE 

In order to further clarify the issue, in this section 
we construct another dynamical black hole solution with 



7 



phantom scalar field. To this end, let us first recall the 
Fonarev solution 13511 which is derived from the action 



S = J d 4 x^g 



R - \d^<j> - V (<j>) 



(41) 



This action describes an Einstein-scalar system. The po- 
tential is assumed to be exponential 



V (0) = V e 



(42) 



describing an Einstein-phantom system with an invariant 
potential. The corresponding Einstein and Klein-Gordon 
equations are 

Gfj, u = -d^d v <f) + -g^ v d a <j>d a (f) - g^V, (49) 



d a d a (j) + ^- = 0. 



(50) 



where Vo and A are two real constants: A describes the 
slope of the potential and Vo = V (0 = 0). Without loss 
of generality, we assume that they are both positive. This 
potential has been investigated at length in the context 
of homogenous cosmology [3^ . The Einstein equation 
and Klein- Gordon equations are 

G^u = d^d v (j) - -gp V d a (pd a <f) - g^V , (43) 



V a V Q c 



dV 



. 



(44) 



The spherically symmetric solution obtained by Fonarev 
is 



ds 2 



= A ln 



-f (r) dr, 2 + 



Vo(\ 2 -2) 2 
2a 2 (6- A 2 ) 



1 



f 2 (r) 



dr 2 + S (r) 2 dn 2 



A In a 



s/X 



2 + 2 V r 



(45) 



where 



Therefore, we construct a dynamical solution represent- 
ing a black hole immersed in a phantom scalar field fluid, 
and corresponding to 



ds 2 



-f (r) 2 dr] 2 + ^r^dr 2 + S (r) 2 dtf 



firy 



1 Vq(A 2 + 2)' 
A 11 2 (A 2 + 6) 
1 



A In a 



where 



lull-^ 

VA^2 V r 



f= 11-2= 

r 



(51) 



A 



S 2 



r 2 I 1 



2w 
r 



VA^2 
1 



(52) 



We assume that A > \J~2~. Next, let us investigate the 
physical meaning of the constant w. When A >> ^/2, we 
have a — I and a = — 1, hence the metric becomes 



/ 
S 2 



2w 
r 



A 



VA^+2' 

2 

i = a |?7| A ' 2 - 2 . 



(46) 



where w and ao are integration constants. The physi- 
cal meaning of w will be clear later, while ao is a nor- 
malization constant for the scale factor. For simplicity, 
hereafter we set ao — 1. When w = 0, the metric (|4"5"]) 
reduces to the FRW universe while, when a = l,a = 1, 
it reduces to the Schwarzschild solution. So this solution 
describes a black hole embedded in FRW universe. 

Now, a dynamical phantom black hole solution can be 
constructed from the Fonarev solution via the transfor- 
mations 



4> — > i(j>, A 
Then, the action becomes 



iX. 



S = d i xy/—g 



(47) 



(48) 



ds 2 



1-^] drj 2 
r 



2 <i- 2 -^) dn 2 



1 - dr 2 
r 



The coordinate transformation 



y = r[l 



2w 



transforms Eq. (f53|) into 



ds 2 



2w 



1 + — - ) <V + ( 1 + — ) dy 



+y 2 dfl 2 , 



(53) 



(54) 



2w 



(55) 



which is the Schwarzschild solution with the mass M — 
—w. Therefore, the parameter w corresponds to the neg- 
ative mass. Hereafter, we use the symbol — M instead of 

UK 
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In order to study whether the singularity can ever be 
naked during the evolution of the phantom field, we cal- 
culate the black hole AH. Using the parameters M and 
a, the phantom black hole solution can be written as 



ds 2 



f 




2a 2 -2 








2J f) 




+r 2 (l 


2M 

r 



2M_y 

dr 2 



d-q 2 



l+a 



dn< 



(56) 



By replacing the conformal time rj with the comoving 
time t, we can rewrite the metric as 



ds 2 



2M_y 



2M 
r 



dt 2 + a (ty 



2M 
r 



dr 2 



a{t) = (to-*)" 



(57) 



The scale factor tells us that the constant to nas the 
physical meaning of the moment at which the Big Rip 
happens. We have a < — 1 because of A > \pl. The 
exponent a is determined by the slope of the potential. It 
is apparent that, when M — 0, the metric (|5T[) describes 
a phantom-dominated Universe. For simplicity, we set 
a = —3, i.e., A = 3/2. Then, the metric (|57[) reduces to 



ds 2 



2M 
r 



-3 



dt 2 



a(ty 



2M 



+r 2 1 



dW 



a(t) = (t -t)- 16/9 . 

Using the physical radius R = ar(l 
of thccoordinate radius r, we obtain 



(58) 

2M/r)~ 1 instead 




(59) 



as the equation of the AH. Here, H = is the Hubble 
parameter. By setting x = 1 + •> /l 



dt 



one obtains 



aM Hx - Ax 2 + 12x - 8 = 



(60) 




-2.185 -2.18 -2.175 - 
t 



2.165 



FIG. 2: The black hole apparent horizon (solid line) increases, 
while the cosmic apparent horizon (dotted line) decreases with 
the divergence of the phantom field. There is a moment at 
which the two horizons coincide and after which they disap- 
pear, the black hole singularity becoming naked. These plots 
correspond to a = —3, to = 1, M = 1. 



We find that Eq. ([60]) has only two physical roots which 
denote the cosmic AH R c and the black hole AH Rf,, re- 
spectively. Using Eq. (|60[) and the definition of x, we 
have plotted the relation between the comoving time t 
and the AHs in Fig. (2). The latter shows that the black 
hole apparent horizon increases while the cosmic appar- 
ent horizon decreases with the blow-up of the phantom 
scalar field. There is a moment at which the two horizons 
coincide and after which both horizons disappear and the 
black hole singularity becomes naked. 



VI. DISCUSSION AND CONCLUSIONS 

The well-known McVittie solution can not describe a 
black hole embedded in a FRW universe. On the con- 
trary, the new solutions of [l7| do: the surface r = Mq/2 
is not an event horizon but a 2-sphere contained inside 
the black hole AH and a conformal Killing horizon when 
M(t) = Moa(t). The energy density and pressure are 
regular at the AH, contrary to the case of the McVittie 
metric in which the pressure diverges on the would-be 
horizon. Therefore, we adopt the solution p7|) of Faraoni 
and Jacques [13 to investigate the problem of phantom 
matter accretion by black holes. We take into account the 
backreaction, which was ignored in the previous work of 
Babichev et al. In a low density background this 

effect can be safely neglected; however, if the energy 
density of the background is large {e.g., if it becomes 
comparable to the black hole density 1/Mq), the metric 
describing this black hole will be modified significantly. 
Therefore, backreaction must be part of the picture. 

Using this solution we find that, according to both the 
quasi-local and the generalized Misner-Sharp mass no- 
tions, the black hole mass increases during the cosmic 
evolution with accretion of the phantom fluid. First of 
all, one notes that more than one notion of mass can be 
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employed, and it is not trivial to identify the correct one 
(the Hawking-Hayward quasi-local mass has proved to be 
the most convenient in previous studies this is- 

sue is absent when using a constant mass Schwarzschild 
black hole as in [8J. The discrepancy between our re- 
sult of increasing Hawking mass and size of the AH, and 
the decreasing mass property found in [f| might be ex- 
plained by the fact that the quasi-local mass does not 
receive contributions from the negative pressure. 

In Faraoni and Jacques [13], there is also a solution 
(eq. (98)) which yields the opposite result, i.e., the black 
hole mass decreases with the accretion of phantom fluid. 
However, that solution requires an outgoing radial en- 
ergy flux q — — (P + p)u/2, which is not easy to interpret 
for an accreting black hole, and could easily explain the 
decreasing mass. Here, we have shown that, for an in- 
ward energy flux, the black hole mass increases even when 
phantom energy fluid is accreted. 

Unfortunately, it is not easy to compare our solution 
to the Babichev et al. solution 0] directly, as they have 
ignored the backreaction effect of expansion, while we 
have simplified our model of fluid flow instead. The 
comparison is not straightforward: it is not a simple 
matter of taking a limit in the relevant equations, be- 
cause here the relevant black hole mass is taken as the 
Hawking-Hayward quasi-local mass for the cosmological 
black hole immersed in a fluid that sources the surround- 
ing FRW universe-this is forced upon us when looking 
for a meaningful mass concept to use. When compar- 
ing with [8j, one must replace this quasi-local mass with 
the Schwarzschild mass and neglect its time dependence, 
which implies setting a to a constant and, in turn, means 
that the background fluid becomes a test fluid (it no 
longer causes the cosmic expansion). The crucial dif- 
ference between our solutions and that of 0] is that the 
black hole mass is no longer obtained as a "mass con- 
trast" with respect to the background, but it will include 
both energy density p and pressure P inside the hori- 
zon, thus causing a flow of accreted total energy that 
is negative (however, in this regime, the effect on the 
Schwarzschild mass is neglected) . While we do not doubt 
that the analysis of @ is correct in the test-fluid regime, 
it is its extrapolation to a gravitating fluid and to an 
expanding universe that we question, and we can only 
propose some special classes of exact solutions to further 
the analysis. Because these may be very special, we do 
not claim that our is the last word on this issue - it is 
an intermediate step (but a valuable one because exact 
solutions are rare), and we auspicate that more exact 
solutions incorporating an explicit radial mass flow be 
found to discuss this issue further. 

The central, and most interesting point of our study, 
however, is not the behavior of the black hole mass but 
rather the process in which the black hole AH expands 
while the cosmic AH shrinks. When the comoving time 
derivative of the scale factor approaches the critical value 
1/(8ttM ) at the time given by Eq. (P?)) . the two hori- 
zons coincide. At times t > t* both horizons disappear 



and the black hole singularity becomes a naked singu- 
larity; this happens before the universe reaches the Big 
Rip. The shrinking of the cosmic apparent horizon has 
serious conse quen ces. The Bekenstein-Hawking formula 
S = A/A [37l l38l| tells us that the horizon entropy S is 
proportional to the area A of its cosmic apparent horizon. 
So, the entropy of a phantom energy-dominated universe 
decreases, in violation of the second law of thermody- 
namics. If the black hole singularity becomes naked in a 
finite time well before the Big Rip, the Cosmic Censor- 
ship Conjecture is violated. By combining these consid- 
erations, one could go as far as reaching the conclusion 
that nature may not permit the existence of phantom 
energy. There is, however, another consideration to keep 
in mind: the solution (| 1 7[) (as well as the similar solu- 
tions of [It], [H| ) is meant to describe a black hole that 
has existed forever, not one that is created by gravita- 
tional collapse. According to the expression of the 
AH radius on which we base our considerations, if this 
black hole is created together with the universe in a Big 
Bang singularity, early on the derivative a was arbitrar- 
ily large and the central singularity was not covered by 
a horizon, either. Hence, the metric (TT7)) describes a 
naked singularity (violating the Cosmic Censorship Con- 
jecture) which gets dressed by a horizon shortly after the 
Big Bang, and evolves as a black hole spacetime until 
it becomes again a naked singularity (violating Cosmic 
Censorship again). This bizarre behavior leads one to 
wonder whether this is a physically meaningful solution. 
After all, it is well known that naked singularities form 
during simulations of gravitational collapse, but they are 
not generic. "Typical" choices of initial data result in 
the collapse to a black hole and naked singularities ap- 
pear to be the exception rather than the rule. It could 
well be that the solution (TIT)) describes a non-physical, 
fine-tuned situation with bizarre phenomenology. The 
more general solutions considered in Sec. IV exhibit the 
same phenomenology and this fact tends to exclude fine- 
tuning, but no truly definitive statement can be made at 
present to this regard due to the paucity of exact solu- 
tions describing the physical situation at hand. 

Using Eq. ([T3]) and Eq. (J23J) with w > — 1, one finds 
that the black hole AH first shrinks and then expands, 
while the cosmic AH always expands with the evolution of 
the universe, with the cosmic AH expanding faster than 
the black hole AH. There was a moment at which the 
two horizons coincided and before which the singularity 
was naked. As an example, we have plotted in Fig. [3] the 
evolution of the cosmic and black hole AHs in a dust- 
dominated universe. 

The reason for the black hole AH shrinking initially is 
that the density, and thus the expansion speed of the uni- 
verse, were very large in the vicinity of this coincidence. 
In this sense, the expansion of the universe "wins" over 
the attraction of the black hole near this coincidence. 

Finally, we comment on the global structure of the cos- 
mological black holes discussed. In Figures 4 and 5 we 
plot the Penrose-Carter diagrams for both the phantom- 
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FIG. 3: The black hole apparent horizon (lower line) first de- 
creases and then increases, while the cosmic apparent horizon 
(upper line) always increases in the dust-dominated Universe. 
There was a moment at which the two horizons coincided and 
before which the singularity was naked. The plots correspond 
to M = 0.8, a = 1. 
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FIG. 5: The Penrose diagram for a dust-dominated universe 
with a black hole. This diagram also coincides with the one for 
the Schwarzschild solution, except for two sections removed 
and two cosmic singularities added. We only label the black 
hole singularity (BHS), white hole singularity (WHS), Big 
Bang singularity (BBS), and Big Crunch singularity (BCS) 
in the diagram. The black hole apparent horizon (BAH) and 
cosmic apparent horizon (CAH) are also shown. 
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FIG. 4: The Penrose-Carter diagram for a phantom- 
dominated universe with a black hole. The diagram is sim- 
ilar to the one for the Schwarzschild solution except for two 
sections removed and two cosmic singularities added. For 
brevity, we only label the black hole singularity (BHS), white 
hole singularity (WHS), Big Bang singularity (BBS), and Big 
Rip singularity (BRS) in the diagram. The black hole appar- 
ent horizon (BAH) and cosmic apparent horizon (CAH) are 
also shown. 



Schwarzschild anti-de Sitter metric is 
2M 



ds z = 



1 - 



r 

2M 



H 2 r 2 j dt 2 
% -l 

2„2 ' 



H 2 r 



dr 1 



r 2 dn 2 . (61) 



In this case, J? is a constant and the black hole AH is 
smaller than that of the corresponding asymptotically 
flat black hole. Accordingly, the generalized Misner- 
Sharp mass is smaller than the one of the hole without 
cosmological constant. We note that the Schwarzschild- 
anti-de Sitter spacetime has a remarkable difference in 
comparison with Schwarzschild-de Sitter: its energy den- 
sity p — — -^H 2 is negative. This solution suggests that 
if a black hole is surrounded by negative and decreasing 
energy density, its mass will decrease monotonically and 
will tend to zero in a finite or an infinite cosmic time. 
However, even though in the solutions considered here 
the phantom kinetic energy density is negative, its total 
energy density is positive, hence this would not happen. 



dominated and the dust-dominated universe with a black 
hole, respectively. The global structure and Penrose 
diagram of the Sultana-Dyer solution has been stud- 
ied in 39], and these are similar for our conformally 
Schwarzschild solutions in a dust-dominated universe (see 
Fig. 3 of Hi). 

Hence, the singularity will never become naked in a 
universe dominated by ordinary (or dark) matter. In a 
universe driven by exotic phantom matter, the black hole 
AH also expands. One then wonders which mechanism 
can make the black hole AH to shrink. Recall that the 
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